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We  show  that  there  are  shape-independent  upper  bounds  to  the  extinction  cross  section  per  unit  volume 
of  dilute,  randomly  arranged  nanoparticles,  given  only  material  permittivity.  Underlying  the  limits  are 
restrictive  sum  rules  that  constrain  the  distribution  of  quasistatic  eigenvalues.  Surprisingly,  optimally 
designed  spheroids,  with  only  a  single  quasistatic  degree  of  freedom,  reach  the  upper  bounds  for  four 
permittivity  values.  Away  from  these  permittivities,  we  demonstrate  computationally  optimized  structures 
that  surpass  spheroids  and  approach  the  fundamental  limits. 
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Many  applications  [1-8]  employ  disordered  collections 
of  particles  to  absorb  or  scatter  light,  and  the  extinction 
for  a  given  total  particle  volume  (for  a  dilute  system  in 
which  coagulation  and  multiple  scattering  are  negligible)  is 
determined  by  the  total  extinction  (scattering  +  absorption) 
cross  section  per  unit  volume  o^^jV  of  the  individual 
particles  [9,10].  In  this  Letter,  we  prove  fundamental  upper 
bounds  on  <7^x1/ L  for  small  metallic  particles  of  any  shape, 
we  show  that  previous  work  on  maximizing  particle 
scattering  [10-14]  (including  “superscattering”  [15-17]) 
was  a  factor  of  6  or  more  from  these  bounds,  and  we 
employ  a  combination  of  analytical  results  and  large-scale 
optimization  (“inverse”  design)  to  discover  nearly  optimal 
particle  shapes.  Most  previous  work  in  this  area  was 
confined  to  spheres  [10,12,16,18,19]  or  a  few  high- 
symmetry  shapes  [11,13-15,17],  whereas  we  optimize 
numerically  over  shapes  with  ;^1000  free  parameters 
(and  prove  our  theorem  for  completely  arbitrary  shapes) 
over  the  visible  spectrum,  and  we  also  consider  coated 
multimaterial  shapes.  We  find  that  the  optimal  is 

invariably  obtained  for  subwavelength  particles  where 
absorption  dominates  and  the  quasistatic  approximation 
applies.  We  can  then  apply  a  little-known  eigenproblem 
formulation  of  quasistatic  electromagnetism  in  terms  of 
“resonances”  in  the  permittivity  e  {not  in  the  frequency  co) 
[20-24],  and  we  employ  various  sum  rules  of  these 
resonances  [20,25,26]  to  derive  a  bound  on  the  cross 
section.  Surprisingly,  very  different  optimized  shapes  (such 
as  ellipsoids  or  “pinched”  tetrahedra)  exhibit  nearly  iden¬ 
tical  <7ext(^^)  spectra  (greatly  superior  to  nonoptimized 
particles)  once  <7ext  is  averaged  over  incident  angle,  a  result 
we  can  explain  in  terms  of  the  quasistatic  resonances. 
Finally,  we  explain  how  our  bounds  provide  materials 
guidance  in  various  wavelength  regimes,  with  potential 
applications  ranging  from  cancer  therapy  [1-3]  and 


plasmonic  biosensors  [4-6,27]  to  next-generation  solar 
cells  [28]  and  optical  couplers  [29]. 

Some  previous  bounds  on  optical  properties  of  dilute 
particle  suspensions  have  been  derived  [30-34].  Purcell 
derived  a  sum  rule  limiting  the  integral  over  all  frequencies 
of  extinction  by  spheroids  [35].  The  limit  has  been 
extended  to  a  variety  of  materials  and  structures  [30-33], 
but  is  geometry  dependent  and  difficult  to  apply  as  a 
general  rule.  Alternatively,  many  authors  have  bounded  the 
effective  “metamaterial”  permittivity  of  composite  media 
[36-39],  a  related  but  not  identical  problem.  The  methods 
presented  here,  applied  to  the  effective  permittivity  of  a 
lossless  dielectric,  are  able  to  reproduce  the  well-known 
Hashin-Strikman  bounds  [40,41]  of  composite  theory. 

A  single  numerical  optimization  conceptually  demon¬ 
strates  many  key  findings  for  nanoparticle  extinction.  To 
illustrate,  we  design  a  silver  particle,  for  maximum  fre¬ 
quency-  and  angle-averaged  extinction  cross  section  per  unit 
volume,  <7ext/F,  over  a  33  nm  bandwidth  at  center  wave¬ 
length  X  =  437  nm  {Q  =  13).  We  do  not  impose  quasistatic 
conditions  a  priori',  we  employ  the  full  Maxwell  equations. 
Ultimately,  the  optimizations  always  converged  to  very 
small,  essentially  quasistatic  sizes,  and  there  is  reason  to 
believe  that  such  quasistatic  sizes  may  be  globally  optimal 
for  metals.  For  quasistatic  particles  (size  approaches  0),  a/V 
is  a  constant.  As  the  particle  size  increases  and  moves  away 
from  the  quasistatic  limit,  the  number  of  surface  modes 
increases  proportional  to  the  surface  area,  so  o/V  decreases. 
Further  size  increases  reach  the  geometric-optics  limit, 
where  a/V  ^  0  as  V  ^  cx). 

We  employed  a  number  of  techniques  to  make  the 
optimization  tractable.  To  quickly  solve  Maxwell’s  equations, 
we  used  a  free-software  implementation  [42,43]  of  the 
boundary-element  method  [44].  Angle  averaging  is  essen¬ 
tially  free  with  such  a  solver.  In  many  apphcations,  the  figure 
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FIG.  1  (color  online).  Angle-averaged  extinction  cross  section 
per  unit  volume  of  computationally  optimized  Ag  particles, 
designed  for  2o  =  437  nmandAi  =  33  nm.  An  ellipsoid  provides 
almost  twice  the  extinction  of  an  optimally  coated  sphere,  but 
optimizing  over  1000  spherical  harmonics  basis  functions  yields 
only  a  2%  further  improvement,  due  to  fundamental  limits  on 
the  eigenvalue  distribution.  Surface  coloring  depicts  the  charge 
density  on  resonance,  where  6'Ag(437  nm)  »  — 5.5  +  0.7/.  Particle 
dimensions  are  10  nm. 

of  merit  is  a  frequency-averaged  extinction,  defined  by  the 
integral  (Text  =  /  We  efficiently  compute 

this  integral  by  contour  integration,  which  for  a  Lorentzian  H 
of  bandwidth  Aco  reduces  to  a  single  scattering  problem  at  a 
complex  frequency  odq  +  iAco  [45^7].  Eor  optimization,  the 
particle  shape  was  parametrized  by  the  zero  level  set  of  a  sum 
of  spherical  harmonics  [48],  i.e.,  r{6,  (j))  —  J2im^imYim  {o.4>) 
(restricting  us  to  “star-shaped”  structures).  Given  the  gradient 
of  the  objective  with  respect  to  these  ;^1000  degrees 
of  freedom  (efficiently  computed  by  an  adjoint  method 
[49-51]),  we  employ  a  free-software  implementation  [52] 
of  standard  nonlinear  optimization  algorithms  [53]  to  find  a 
local  optimum  from  a  given  starting  point.  We  also  optimized 
the  few  degrees  of  freedom  of  coated  spheres  and  ellipsoids 
for  the  sake  of  comparison. 

Eigure  1  depicts  the  optimal  particles  and  their  respective 
extinction  spectra.  The  optimal  designs  were  in  the  quasi¬ 
static  limit,  with  dimensions  ;^10  nm.  A  10  nm  size  is  not 
uniquely  optimal,  but  is  rather  the  size  at  which  performance 
is  dominated  by  quasistatic  response.  Doubling  the  size  is 
clearly  worse  (;^2%),  whereas  a  tenfold  size  reduction  is 
better  by  only  0.2%,  within  the  meshing  error.  Eurthermore, 
our  calculations  do  not  account  for  the  quantum  plasmonic 
properties  of  very  small  particles,  which  would  tend  to 
decrease  performance  below  10  nm  dimensions  [54]. 

We  see  that  uncoated  ellipsoids  provide  significant  gains 
over  coated  spheres,  which  already  provide  a  substantial 
response  [1,3,10,55]  (coated  ellipsoids  showed  no  further 
benefit).  This  suggests  a  principle  that  tuning  resonances  by 
geometrical  deformation  rather  than  by  coatings  enhances 
performance.  Oblate  (“pancake”)  ellipsoids  are  superior  to 
prolate  (“rod”)  elhpsoids,  because  they  couple  to  two  of  the 
three  polarizations  of  randomly  oriented  incident  waves.  In 
the  much  larger  spherical-harmonics  design  space,  the  optimal 


structure  turned  out  to  be  a  “pinched”  tetrahedron  (FT),  which 
can  be  conceptualized  as  pinching  a  sphere  towards  the  four 
centroids  of  the  faces  of  an  inscribed  tetrahedron.  Surprisingly, 
the  much  larger  design  space  yielded  a  structure  that  was  only 
2%  better  than  the  best  ellipsoid.  The  two  structures  have  very 
different  responses  for  a  given  incidence  angle  and  polariza¬ 
tion;  only  when  averaged  over  angle  and  polarization  do  the 
responses  become  nearly  identical.  Also  shown  in  Eig.  1  are 
the  imaginary  parts  of  the  charge  densities  for  resonant 
incident  waves,  explained  below.  Intuitively,  the  ellipsoid 
and  PT  are  better  than  a  coated  sphere  because  the  opposing 
surface  charges  have  larger  spatial  separations. 

The  nearly  identical  spectra  for  the  spheroid  and  PT  can 
be  explained  by  a  fundamental  restriction  on  quasistatic 
eigenmodes,  which  are  prevented  from  fully  coupling  to 
external  radiation.  In  the  quasistatic  limit,  the  incident  field  is 
locally  constant  and  the  response  of  the  system  is  determined 
by  induced  charge  densities  at  the  surfaces.  One  can 
constmct  the  fields  from  the  homogeneous  Green’s  functions 
of  the  induced  surface  charges  c-(x).  Eor  a  surface  S,  the 
surface  integral  equation  for  the  charge  density  is  [20-24]: 

A(t{x)- jn{x)-G^{x-x')a{x')dS'  =  E^^{x)-n{x),  (1) 

V* 

ko 

where  A  =  (^j^t  +  ^ext)/2(Gnt  “  ^ext)  relates  interior  and 
exterior  permittivities,  the  electrostatic  Green’s  function 
G^(x)  =  x/4.;r|xp,  and  E^’^^(x)  •  n(x)  is  the  normal  com¬ 
ponent  of  the  incident  field  at  x  (boldface  indicates  vector 
quantities).  As  distinguished  from  the  VQ^ondinifrequencies  of 
Maxwell’s  equations,  there  are  resonant  permittivities  Gnt/^ext 
for  the  quasistatic  integral  equation.  These  are  negative,  real¬ 
valued  permittivities  at  which  self-sustaining  charge 
densities  exist  without  external  fields,  for  specific  eigenmodes 
(T„  satisfying  where  K  is  the  Neumann-Poincare 

integral  operator  defined  by  Eq.  (1).  The  eigenvalues  lie  in 
the  interval  [—  1  /2, 1  /2]  [23,24,56],  such  that  <  0.  The  left 
eigenvectors  of  K,  denoted  have  the  same  eigenvalue 
spectrum  as  the  (T„  (i.e.,  =  /1„t„)  and  provide  the 

orthogonality  condition  =  J^o^TrudS  =  [23]. 

Equation  (1)  is  valid  for  linear,  isotropic,  and  non¬ 
magnetic  materials.  Its  generalization  to  multiple  surfaces 
takes  A  to  a  diagonal  matrix  [57];  the  eigenmode  decom¬ 
position  of  k  imposes  strict  requirements  on  the  allowable 
form  of  the  matrix,  such  that  each  interface  must  separate 
the  same  materials.  Thus,  Eq.  (1)  is  valid  for  arbitrarily 
many  interacting  objects,  possibly  coated  or  holey  (e.g., 
torii),  as  long  as  there  are  only  two  permittivities. 

The  eigenmodes  of  K  contribute  to  absorption  and 
scattering  through  a,  the  particle’s  polarizability  per  unit 
volume  V,  which  relates  the  incident  field  to  the  dipole 
moment  by  p^  =  Decomposing  the  charge 

density  as  a  superposition  of  eigenmodes,  a  =  X]„c„(T„, 
solving  for  via  Eq.  (1)  and  for  the  dipole  moment  via 
p  =  fs  xodA,  yields 
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L„  -  ^{(o)  ’ 


(2) 


where  n„)/V  is  the  dipole  strength  of 

each  mode,  =  1/2  —  is  the  depolarization  factor,  and 
^{co)  =  — qnt/(^int  “  ^ext)  ^presents  the  relative  properties 
of  the  interior  and  exterior  materials. 

The  distribution  of  eigenmodes,  and  therefore  the 
induced  susceptibility,  is  restricted  by  two  crucial  sum 
rules.  The  first  is  the  /-sum  rule  [25,58],  limiting  the  total 
dipole  strength  for  uncoated  particles: 


The  total  dipole  strength  of  coated  particles  is  reduced  by  the 
metallic  volume  fraction  /  [58].  The  second  sum  rule  [26,58], 
apphcable  for  coated  and  uncoated  particles,  states  that  the 
weighted  average  of  the  depolarization  factors  must  be  1/3: 


(Ln) 


'YlnPn^n  _  ^ 

J2nPn  3’ 


(4) 


where  denotes 

A  sphere  has  a  depolarization  factor  of  1/3,  leading  to  a 
‘‘plasmon”  resonance  at  (^  =  1/3).  Equation  (4) 

dictates  that  the  average  depolarization  factor  of  every 
structure  must  equal  that  of  the  sphere.  Although  it  was 
exploited  for  composites  with  certain  symmetries  [59,60], 
this  general  property  has  not  been  widely  recognized  and  is 
very  important  in  limiting  possible  extinction  rates. 

The  average  extinction  of  randomly  arranged  particles  is 
proportional  to  the  imaginary  part  of  Tra^^  [9],  which  is 
given  by  Eq.  (2): 


^ext 


2jt 


Vim 


Ln  -  ^(C0)_ 


Pn- 


(5) 


A  resonance  occurs  for  =  ^^(u;),  where  r  and  i  sub¬ 
scripts  denote  real  and  imaginary  parts,  respectively.  Eor 
particles  in  vacuum  with  susceptibility  x{^)  — 

^{(d)  =  —\/x{ot)).  Only  metals,  with  6^(09)  <  0,  can 
achieve  0  <  <  1,  and  therefore  exhibit  quasistatic  sur- 

face-plasmon  modes.  The  second  sum  guarantees  that 
(except  in  the  case  =  1/3)  a  particle  cannot  have  all 
of  its  dipole  strength  on  resonance;  there  must  always  be  a 
counterbalancing  dipole  moment  such  that  (L„)  =  1/3. 

Eor  a  given  material  parameter  ^^(u;),  we  can  show  that 
the  optimal  distribution  of  eigenmodes  has  at  most  two 
distinct  depolarization  factors,  and  L2.  We  have  rigor¬ 
ously  derived  the  exact  locations  of  the  two  eigenvalues 
[58],  but  for  relevant  materials  a  simple  solution  suffices: 


0<^r<  1/3 

(OXr) 

1/3  <  <  1 

(6) 

(0,1) 

<  0  or  >  1, 

which  corresponds  to  placing  as  much  of  the  dipole  moment 
as  possible  on  resonance  (L  =  ^^),  and  the  rest  of  the  dipole 
strength  at  the  opposite  boundary  to  satisfy  the  second  sum 
rule.  Equation  (6)  is  exact  for  =  0  (both  a  low-loss  =  0 
and  infinite-loss  Xi  00  limit),  but  is  also  very  accurate 
(error  <  10“^)  otherwise.  With  given  by  Eq.  (6),  we  can 
solve  for  from  Eqs.  (3)  and  (4).  Plugging  and  into 
Eq.  (5)  yields  the  upper  limit  to  the  extinction  per  unit  volume: 

2z?(l+Zr)+z/3+2/,+4/5)+2j4  ^  Xr  ^  \ 

Xi{x]H^+Xrf)  IzP  3 

(7) 

which  provides  a  limit  for  any  possible  susceptibihty, 
independent  of  geometry.  Ideal  scatterers  are  uncoated,  and 
have  metallic  permittivities  with  small  imaginary  parts  (as  in 
Ref.  [19])  and  very  negative  real  parts;  for  61  ^  Eq.  (7) 
simplifies  to 


^ext  ^  ^ 


(Text  4;r  6:^  _ ,  , 

-^< - -kO(ei), 


(8) 


where  the  “O”  notation  indicates  the  asymptotic  scaling  of  the 
higher-order  term. 

Equations  (7)  and  (8)  represent  fundamental  limits  to 
quasistatic  particle  extinction.  Eigure  2  illustrates  these  limits 
by  normalizing  them  relative  to  the  value  of  extinction  on 
resonance,  (Tj-es  =  2;i/3A^i{co),  and  comparing  them  to  elhp- 
soid  limits  computed  through  nonlinear  optimization  [52]. 
The  structural  eigenmodes  were  computed  with  boundary- 
element  method  software  [61].  (Text/<^res  can  be  thought  of  as 
the  number  of  fully  coupled  polarizations;  only  at  =  1/3 
(ey.  —2)  can  full  couphng  to  all  three  polarizations  occur. 
Thus  we  see  why  ellipsoids  perform  very  well,  and  why  the 
optimal  structure  of  Eig.  1  barely  outperformed  the  ideal 
elhpsoid:  in  many  cases,  full  couphng  to  two  polarizations 
closely  approaches  the  ideal  performance.  This  is  exactly  true 
for  €y.  —00,  one  of  the  cases  in  which  elhpsoids  reach  the 

upper  bound.  The  other  three  cases  are  =  —2,  =  —1, 

and  ey.  =  0,  for  which  a  sphere,  infinite  cylinder,  and 
infinitely  thin  disk  are  optimal,  respectively.  In  each  case, 
the  spheroid  depolarization  factors  [9]  are  identical  to  those  of 
the  optimal  general  shape,  given  by  Eq.  (6). 

Included  in  Eig.  2  are  optimizations  at  other  permittiv¬ 
ities  (assuming  the  complex  permittivity  of  Ag);  there  is  a 
family  of  “pinched  tetrahedron”  structures  that  emerge  as 
superior  design  choices  over  ellipsoids.  It  is  important  to 
note  that  spheres  are  not  globally  optimal,  as  the  normali¬ 
zation  factor  (Tj-es  is  a  function  of  ey..  The  inset  shows  the 
absolute  extinction,  which  scales  as  el/ei. 

The  limits  of  Eqs.  (7)  and  (8)  may  appear  to  contradict 
arguments  in  coupled-mode  theory  (CMT)  [16,62],  but  in 
fact  do  not.  CMT  predicts  2^  scaling  only  when 

radiation  loss  dominates  over  absorption  loss;  when 


123903-3 


PRL  112,  123903  (2014) 


PHYSICAL  REVIEW  LETTERS 


week  ending 
28  MARCH  2014 


FIG.  2  (color  online).  Fundamental  extinction  limits,  given  by 
Eq.  (7)  and  normalized  to  the  maximum  extinction  of  a  single¬ 
polarization  resonance  Spheres  are  not  optimal  for  absolute 
(Text  (see  inset),  but  do  enable  full  coupling  to  three  polarizations, 
given  by  the  normalized  value  Markers  indicate 

computationally  optimized  structures,  ^/(m)  is  taken  to  be  that 
of  Ag,  although  this  has  only  a  small  effect  on  the  line  shape. 
Ellipsoids  can  approach  the  general  bounds  in  four  limits: 
€r  —00  (oblate  disk),  =  —2  (sphere),  e^.  —  —I  (cylinder), 
and  €r  —  0  (oblate  disk).  Computationally  optimized  pinched 
tetrahedra  improve  upon  ellipsoids  at  intermediate  Inset: 

upper  bound  on  which  increases  with  ej/ei  (e^  <  0). 
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FIG.  3  (color  online).  Distributions  of  dipole  strength  for 
computationally  optimal  [pinched  tetrahedron  (PT)],  nearly 
optimal  (ellipsoid),  and  nonoptimal  (torus  and  bow-tie  antenna) 
structures  subject  to  the  bounds  of  Eq.  (7).  The  PT  and  ellipsoid 
have  degenerate  modes  at  the  optimization  permittivity,  in  this 
case  €r{o))  =  —3.2.  The  PT  outperforms  the  ellipsoid  because  its 
“undesirable”  eigenmodes  are  closer  to  L  =  1,  enabling  a  larger 
dipole  strength  at  Lj  =  0.24  (ej  =  —3.2).  Equally  important  is 
the  lack  of  other  bright  modes;  e.g.,  torii  (crosses)  and  bow-tie 
antennas  (open  circles)  have  disperse  modes,  reducing  overall 
extinction.  The  PT/ellipsoid  modes  coincide  at  Lj  =  0.24  but  are 
split  for  visualization. 


absorption  dominates,  CMT  predicts  as  in 

Eqs.  (7)  and  (8).  Absorption  loss  dominates  for 
(InajXf  ^  El  (Refs.  [18,63]),  which  is  satisfied  by  all 
quasistatic  metallic  particles  in  the  visible  and  infrared. 

Figure  3  shows  the  depolarization  factor  distributions  of 
the  ideal  pinched  tetrahedron  and  ellipsoid  structures,  as 
well  as  nonideal  structures.  We  see  that  the  dipole  moments 
are  largely  concentrated  at  the  desired  permittivity,  except 
as  required  to  keep  the  centroid  of  equal  to  1/3.  The 
tetrahedra  have  the  off-resonance  dipole  moments  distrib¬ 
uted  closer  to  the  boundary  =  1  than  ellipsoids,  explain¬ 
ing  the  slightly  superior  performance. 

Figure  4  illustrates  the  general  utility  of  the  bounds  of 
Eq.  (7).  For  a  given  permittivity,  a  maximum  extinction  per 
unit  volume  can  be  computed  independent  of  structure.  This 
has  important  imphcations  for  material  selection,  which 
varies  by  apphcation  and  frequency.  Although  the  bounds 
are  quasistatic,  as  discussed  earlier  the  quasistatic  bound  may 
be  globally  optimal.  Indeed,  the  infrared  extinction  limits  are 
3  orders  of  magnitude  larger  than  the  best  nonquasistatic 
particles  investigated  to  date  [7].  Although  the  bounds  are  for 
a  single  frequency,  through  complex-frequency  calculations, 
or  known  material  quality  factors  (geometry-independent 
[64]),  rational  design  for  any  bandwidth  is  possible. 

We  can  compare  our  structures  to  recently  proposed  “super¬ 
scattering”  structures  [15-17].  Of  primary  importance  is  the 
figure  of  merit  (FOM).  For  apphcations,  volume  or  weight  is 
the  relevant  normalization.  Normalizing  by  2^,  as  in  [15-17], 
favors  larger  particles  approaching  wavelength  scale.  A 
smaller  particle  with  larger  likely  cannot  extinguish 

a  full  square  wavelength.  Yet  a  dilute  mixture  of  such  particles 


could,  with  much  smaller  volumes.  As  an  example,  tw;o 
quasistatic  nanoellipsoids,  with  an  8:1  major  to  minor  axis 
ratio  can  achieve  the  same  as  the  single  particle  in  [  1 6] , 

while  requiring  1  / 270th  of  the  volume.  A  single  “channel”  in  a 
nonspherical  structure  can  extinguish  much  more  strongly  than 
multiple  channels  in  a  spherical  structure. 

Small,  absorbing  nanoparticles  show  promise  for  a  variety 
of  scientific  and  technical  apphcations.  Experimentally 
approaching  the  limits  derived  here  would  already  represent 
a  significant  achievement.  A  possible  further  improvement 
could  come  from  harnessing  exotic  material  systems  [55,66], 


EIG.  4  (color  online).  Shape-independent  fundamental  limit  to 
extinction  per  unit  volume  for  the  highest-performing  metals  [65] 
at  visible  and  infrared  wavelengths.  A  mixture  of  A1  and  Ag 
nanoparticles,  properly  designed,  could  provide  ideal  extinction 
over  the  visible  and  near-  to  midinfrared.  Inset:  minimum  volume 
fraction  V j)?  required  for  (jgxt  =  2^.  It  is  possible  to  achieve  2^ 
cross  sections  for  V j)?  <  10“^. 
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where  geometry-dependent  material  resonances  cannot  be 
modeled  with  bulk  permittivities. 
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